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We report a search for the decays B~ — > K + tt~ti~ and B~ — * K~ K~n + , which are highly sup- 
pressed in the standard model. Using a sample of (467 ± 5) x 10 6 BB pairs collected with the 
BABAR detector, we do not see any evidence of these decays and determine 90% confidence level 
upper limits of B(B~ -f K + ir~ir~) < 9.5 x 10~ 7 and B(B~ -»■ K~K~i: + ) < 1.6 x 10~ 7 on the 
corresponding branching fractions, including systematic uncertainties. 
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The decays B~ — > K + ti~ii~ and B~ — » _ftT~i ; 4r~7r + 
proceed via 6 — > <icfs and 6 — > ssd quark transitions, re- 
spectively. These are highly suppressed in the standard 
model (SM). Compared with the penguin (loop) transi- 
tions b —>■ qqd, and b — ► qqs they are additionally sup- 
pressed by the small Cabibbo-Kobayashi-Maskawa ma- 
trix 0, % element factor \V td V* s \ ~ 3 x 10~ 4 , leading 
to predicted branching fractions of only O(10" 14 ) and 
O(10~ n ), respectively 0, Q. Example SM decay dia- 
grams can be seen in Fig. [T] 



u,c,t 




W < W 




K + 




FIG. 1: Example standard model decay diagrams for the 
decays B~ —> K + tt~-k~ and B~ — > K~ K~ tt + , respectively. 

These branching fractions can be significantly en- 
hanced in SM extensions such as the minimal supersym- 
metric standard model (MSSM) with or without con- 
served R parity, or in models containing extra U(l) gauge 
bosons. For example, the branching fraction for the 
b — > dds transition B~ — > ir~K*° can be enhanced from 
about 1CT 16 in the SM to about 1CT 6 in the presence of 
an extra Z' boson Q. The branching fraction for the 
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b — > ssd decay B — > K K*° can be enhanced from 
about 7 x 10" 14 in the SM to about 6 x 1CT 9 in the 
MSSM [|. 

Observations of the decays B~ — > K + tt~it~ and 
B~ — > K~ K~n + would be clear experimental signals for 
the b — > dds and b — ► ssd quark transitions 0, [? L These 
decay modes have been previously searched for §|, 
Hlj ]. and the most restrictive 90% confidence level exper- 
imental upper limits B(B~ — ► K + tt~it~) < 1.8 x 10~ 6 
and B(B- -> K-R--K+) < 1.3 x 10~ 6 [l]| were obtained 
from an analysis of 81.8 fb _1 of BaBAR data. Upper lim- 
its on b — > ssd transitions have also been set using the 
decays B~ _— » K*-K~tt + [3, B° -> K*°K-tt + [3, 
and B° -> K*°K*° Q. 

We report the results of a search for the decays 
B~ — > _ftT + 7r~7r~ and B~ — > K~-KT~7r + . Inclusion of the 
charge conjugate modes is implied throughout this paper. 
The data used in this analysis, collected at the PEP-II 
asymmetric energy e + e~ collider [15| . consist of an inte- 
grated luminosity of 426 fb -1 recorded at the T(45) reso- 
nance. In addition, 44 fb _1 of data were collected 40 MeV 
below the resonance and are used for background charac- 
terization. These samples are referred to as on-resonance 
and off-resonance data, respectively. The on-resonance 
data sample contains (467 ± 5) x 10 6 BB pairs. Beyond 
the larger sample size, we utilize improved analysis tech- 
niques for background rejection and signal identification 
compared with our previous study [Iol |. 

The BaBAR detector is described in detail else- 
where [l6[ . Charged particles are detected and their mo- 
menta measured with a five-layer silicon vertex tracker 
and a 40-layer drift chamber (DCH) inside a 1.5 T 
solenoidal magnet. Surrounding the DCH is a detec- 
tor of internally reflected Cherenkov radiation. Energy 
deposited by electrons and photons is measured by a 
CsI(Tl) crystal electromagnetic calorimeter. 

We select B~ — > K + ir~ir~ candidates by combining 
a charged kaon candidate with two charged pion candi- 
dates, each of which has charge opposite to the kaon. 
Similarly, B~ — * K~K~ir + candidates are selected by 
combining two charged kaon candidates with a charged 
pion candidate. Each track is required to have a mini- 
mum transverse momentum of 50 MeV/c, and to be con- 
sistent with having originated from the interaction re- 
gion. Identification of charged pions and kaons is accom- 
plished using energy loss (dE/dx) information from the 
silicon vertex tracker and DCH, and the Cherenkov angle 
and number of photons measured in the detector of in- 
ternally reflected Cherenkov radiation. The efficiency for 
kaon selection is approximately 80% including geometri- 
cal acceptance, while the probability of misidentification 
of pions as kaons is below 5%. The corresponding effi- 
ciency and kaon misindentification rate for the pion selec- 
tion criteria are 95% and less than 10%, respectively. We 
require all charged particle candidates to be inconsistent 
with the electron hypothesis, based on a cut-based selec- 
tion algorithm that uses information from dE/dx, shower 
shapes in the electromagnetic calorimeter, and the ratio 
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of the shower energy and track momentum. 

To avoid a potentially large source of background aris- 
ing from B decays mediated by the favored b — ► c tran- 
sition, we veto B candidates for which pairs of daughter 
tracks have invariant mass combinations in the ranges 
1.76 < m Kn < 1.94GeV/c 2 , 2.85 < m K7T < 3.25GeV/c 2 , 
and 3.65 < ttik-k < 3.75GeV/c 2 . These remove events 
containing the decays D° — ► K~ir + , J/ip — > £ + £~, and 
i/j(2S) — > £ + £~ , respectively, where the leptons in the J/0 
and V , (2S') decays are misidentified as pions or kaons. 

Continuum e + e~ qq (q = u, d, s, ) events are the 
dominant background. To discriminate this type of event 
from signal, we use a neural network that combines five 
variables. The first of these is the ratio of L-i to Lq, with 
Lj = P* \ cos 8* \ 3 , where p* is the particle momen- 
tum, 6* is the angle between the particle and the thrust 
axis determined from the B candidate decay products, 
the sum is over all tracks and neutral clusters not asso- 
ciated with the B candidate, and all quantities are cal- 
culated in the e + e~ center-of-mass (CM) frame. The 
other four variables are the absolute value of the co- 
sine of the angle between the B direction and the beam 
axis; the magnitude of the cosine of the angle between 
the B thrust axis and the beam axis; the product of 
the B candidate's charge and the output of a multivari- 
ate algorithm that identifies the flavor of the recoiling 
B meson [ljj ; and the proper time difference between 
the decays of the two B mesons divided by its uncer- 
tainty. The angles with respect to the beam axis are 
calculated in the CM frame. The neural network out- 
put NN out is distributed such that it peaks around for 
continuum background and around 1 for signal. We re- 
quire NN out > 0.5 (NNout > 0.4) for B~ -> K + -k~-k~ 
(B~ — > K~ K~ir + ) candidates. These requirements re- 
tain approximately 90% of the signal, while rejecting ap- 
proximately 80% of the continuum background. 

In addition to the neural network output, we distin- 
guish signal from background events using two kinematic 
variables: the difference AE between the CM energy of 
the B candidate and y/s/2, and the beam-energy substi- 
tuted mass ities = -\/s/4 — p^ 2 , where y/s is the total 
CM energy and is the momentum of the candidate 
B meson in the CM frame. The AE distribution peaks 
near zero with a resolution of around 19 MeV, while the 
itjes distribution for signal events peaks near the B mass 
with a resolution of around 2.4MeV/c 2 . We select signal 
candidates that satisfy 5.260 < toes < 5.286 GeV/c 2 and 
\AE\ < 0.100 GeV. This region includes a sufficiently 
large range of toes below the signal peak to determine 
properties of the continuum distribution. 

The efficiency for signal events to pass the se- 
lection criteria is 21.6% (17.8%) for B~ -> K^tt~tt~ 
(B~ — > K~ K~n + ), determined with a Monte Carlo 
(MC) simulation in which the decays are generated 
uniformly in three-body phase space. The BABAR 
detector Monte Carlo simulation is based on GEANT 
4 [3 and EvtGen [l|. We find that 8.2% (5.1%) of 
B~ — > K + ir~ir~ (B~ —> K~ K~tt + ) selected events con- 



tain more than one candidate, in which case we choose 
the one with the highest probability that the three tracks 
originate from a common vertex. 

We study possible residual backgrounds from BB 
events using MC event samples. Backgrounds arise from 
decays with topologies similar to the signal but with 
some misreconstruction. Such effects include kaon/pion 
misidentification, the loss of a soft neutral particle, and 
the association of a particle from the decay of the other B 
in the event with the signal candidate or vice versa. We 
find that the backgrounds can be conveniently divided 
into five categories for both the K + ir~ir~ and K~K~w + 
channels, each of which is dominated by one or two par- 
ticular decays but also includes other decay modes that 
result in similar toes and AE shapes. Table U provides 
details of the composition of the background categories. 

In order to obtain the B~ — > K + ir~ir~ and 
B~ — > K~ K~ir + signal yields, we perform unbinncd 
extended maximum likelihood fits to the candidate 
events using three variables: toes, AE, and NN out . For 
each event hypothesis j (signal, continuum background, 
or one of the five BB background categories), we define 
a probability density function (PDF) 



W = P i (m B s i ,A£ i ).p J -(NrC lt ) : 



(1) 



where i denotes the event index. For the signal, con- 
tinuum background, and the BB background categories 
with small correlations between to-es and AE, the PDF 
is further factorized 

Vj(rriES i ,AE i ) = Potties') ■ Pj {AE 1 ) . (2) 
The extended likelihood function is 



C = cxp (^-Y^n^j Y[ 



E 



(3) 



where nj (njt) is the yield belonging to the event hypoth- 
esis j (k). 

The signal toes and AE shapes are parametrized 
with the sum of a Gaussian and a Crystal Ball func- 
tion [2fj, |2lj, l22| and the sum of two Gaussian functions, 
respectively. We determine the shape parameters by tak- 
ing the values obtained from signal MC and correcting for 
differences between data and MC seen in a control sample 
of B~ -> D°tt- with D° -> K-TT+ decays. The contin- 
uu m back ground to-es shape is described by the function 
xyl — x 2 exp [— £[1 — x 2 )], with x = 2toes/-\A and £ a 
free parameter [23| . while the continuum AE shape is 
modeled with a linear function. We describe the toes 
and AE shapes of each BB background category using 
either independent ID histograms or a 2D histogram de- 
termined from MC samples. The decision to use ID or 2D 
histograms is made based on the magnitude of the cor- 
relations between these variables for each category and 
the effect on the signal yield of neglecting such correla- 
tions, discussed below. The PDFs for categories 1, 2, 
and 3, for both B~ -> K+^tt- and B~ -> K-R-TT+, 
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TABLE I: Summary of the B background categories, giving the dominant decay mode, numbers of expected and observed 
events and the character of the jties and AE distributions. "Peaking" indicates that the shape is similar to that of the 
signal. "Broad peak," "left peak," and "right peak" differ from the signal in being wider or shifted to lower or higher values, 
respectively. The number of expected and observed events are also given for the continuum background. 





Category 


1 


2 


3 


4 


5 


Continuum 


Dominant mode(s) 


B~ -> D u tt~; 


B~ — » TT~TT + TT~ 


B~ -> R-TT+TT-&: 




Generic BB 




D° -> K~K+ 




B ^K+7T-7T 








Number of expected events 


80 ±3 


57 ±4 


472 ± 24 


43 ±1 


917 ±19 


25552 ± 495 


Number of observed events 


61 ±70 


-153 ±94 


1116 ±347 


-26 ± 152 


197 ± 273 


25261 ± 198 


mES Structure 


Peaking 


Peaking 


Broad peak 


Broad peak 


Continuum-like 




AE Structure 


Left peak 


Right peak 


Broad peak 


Right peak 


Continuum-like 




B~ -> K-R-TT+ 


Category 


1 


2 


3 


4 


5 


Continuum 


Dominant mode(s) 


B~ -> K~ K + K~ 


B~ -> K~n + n~ 


B~ -> D u n~; 


Generic B + B~ 


Generic B U B U 








D° -> K-it+tt 








Number of expected events 


190 ±9 


198 ±9 


61 ±4 


312 ± 11 


173 ±8 


6088 ± 241 


Number of observed events 


213 ±41 


240 ± 37 


-34 ±55 


380 ± 117 


95 ± 107 


6953 ± 100 


jties Structure 


Peaking 


Peaking 


Broad peak 


Broad peak 


Continuum-like 




AE Structure 


Left peak 


Right peak 


Left peak 


Continuum-like Continuum-like 





arc modeled using 2D histograms. We use ID histograms 
to describe all NN out distributions. These histograms 
are obtained from MC samples for the signal and BB 
background categories, and from a combination of on- 
resonance data, in a continuum-dominated sideband of 
mES and AE, and off-resonance data for the continuum 
background. 

The nine free parameters in our fits are the yields of the 
signal, continuum and all five BB background categories, 
the £ parameter of the continuum mES shape, and the 
slope of the continuum AE shape. 

We test the fitting procedure by applying it to ensem- 
bles of simulated experiments where events are generated 
from the PDF shapes described above for all seven cat- 
egories of events. We repeat the exercise with qq events 
generated from the PDF while signal events are randomly 
extracted from the MC samples. The BB background 
events are either generated from PDF shapes or drawn 
from MC samples. In all cases, these tests confirm that 
our fit performs as expected, with very small biases on 
the fitted signal yields, for which we correct the measured 
yields and include systematic uncertainties. 

We apply the fit described above to the 26 478 
B~ -> K+tt-tt- and 7 822 B~ -> K~K--k + candidate 
events selected from the data recorded at the T(4S) res- 
onance. We find 22 ± 43 and —26 ± 19 signal events, 
respectively, (statistical uncertainties only). The yields 
of continuum and all BB background categories (shown 
in Table [l| are generally consistent with expectations. 
The yields of BB background categories 3 and 5 from 
the fit to B~ — > K + ir~ir~ candidates do not show per- 
fect agreement; however, the sum of their yields is consis- 
tent with the expectation and, owing to the strong neg- 
ative correlation between the yields of these categories, 
the discrepancy with the expectation is not significant. 



Such behavior was seen in the fit validations and has been 
shown not to effect the signal yield. The results of the 
fits are shown in Fig. 

We determine the branching fractions for 
B~ — > K + tt~tt~ and B~ — > K ~ K~tt + by applying 
corrections for the small biases evaluated in the MC 
studies (3.6 ± 2.4 and 0.5 ±1.0 events, respectively) 
and then dividing by the selection efficiencies and the 
total number of BB pairs in the data sample. We 
assume equal decay rates of T(4S) — > B + B~ and 
B°B°. Systematic uncertainties on the fitted yields 
arise from uncertainties in the PDF shapes (10.0 and 
3.5 events, respectively) including possible data/MC 
differences in the signal PDF shapes studied using 
the B~ — > D°ir~ control samples discussed above. We 
estimate the uncertainty on the fit bias (3.1 and 1.0 
events, respectively) to be half the value of the correc- 
tion combined in quadrature with the precision with 
which the bias is known. Uncertainties on the efficiency 
arise from possible data/MC differences for tracking 
(1.2%) and particle identification (4.2%). We consider 
two sources of uncertainty related to the Dalitz plot 
distributions of the signal decays. The first is related 
to the variation of the efficiency over the parts of the 
Dalitz plots that are included in the analysis: from 
MC studies, the uncertainties are found to be 13.0% 
for B~ -> K+ir-n- and 13.5% for B~ -> K-K~tt+. 
The second is due to the correction for the vetoed 
parts of the Dalitz plots, which we estimate for various 
signal decay distributions. In addition to the nominal 
phase-space distribution, we consider decays dominated 
by the intermediate states if* (892) or if,* (1430) 
(modeled using the LASS [24| shape, as implemented 
in our Dalitz plot analysis of B + — > K + ti~it + [25(). 
We mimic a possible enhancement at low n~7r~ or 



5.26 5.27 5.28 „ -0.1 -0.05 0.05 0.1 0.6 0.8 1 
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FIG. 2: Projections of the selected events with the fit results overlaid. The top (bottom) set of plots are for 73" — > K + tv~-k~ 
(B~ — > K~ K~-w )■ From left to right the plots show the projections onto rn.ES, AE, and the output of the neural network. 
The black points are the data, the solid blue curve is the total fit, the dotted red curve is the continuum background, the dashed 
green curve is the total background, and the dash-dotted black curve at the bottom is the signal. The continuum component 
has been suppressed in these plots by applying an additional requirement on the ratio of the signal likelihood to the sum of the 
signal and continuum likelihoods, calculated without use of the plotted variable. The value of the requirement for each plot is 
chosen to reject about 95% of the continuum background while retaining about 55% of the signal. 



K~K~ invariant mass by employing an ad hoc doubly 
charged scalar resonance with mass 1500 MeV/c 2 
and width 300 MeV/c 2 . The efficiency of the veto 
requirement is larger than that for the phase-space 
MC in all alternative models, so we assign asymmetric 
systematic errors of ^"g% for B~ — ► K + Tr~ir~ and 
for B~ — > K~ K~ir + . The uncertainty on the number of 
BB pairs is 1.1%. Including all systematic uncertainties, 
we obtain the following results for the branching frac- 
tions: B(B~ -► #+71-71--) = (1.8 ±4.3 ±0.9) x 10~ 7 and 
B(B- -> K-R--K+) = (-3.2 ± 2.3 x 10~ 7 , where 
the first uncertainties are statistical and the second are 
systematic. 

We have also calculated the branching fractions using 
event-by-event efficiencies ap plied to signal weights ob- 
tained from the fit result [2fJ, |27| • We obtain results con- 
sistent with our main results within the efficiency varia- 
tion systematic uncertainty. We have also checked that 
removing each of the discriminating variables from the 
fit, in turn, gives consistent results. 

To obtain 90% confidence level upper limits on the 
branching fractions, we use the frequentist approach of 
Feldman and Cousins 28]. We determine 90% confi- 
dence region bands that relate the true values of the 
branching fractions to the measured numbers of sig- 
nal events. These bands are constructed using the re- 
sults of MC studies that account for relevant biases 
in the fit procedure and include systematic uncertain- 



ties. The construction of the confidence region bands is 
shown in Fig. [3] The 90% confidence level upper limits 
are found to be B(B~ — > K + ir~n~) < 9.5 x 10~ 7 and 
B(B- -> < 1.6 x 10~ 7 . To aid comparison 

with other experiments, we also extract the sensitivities 
Bq defined as the 90% confidence level upper limits that 
would be obtained in the case of zero fitted signal yield. 
The sensitivities are B (B~ — > A' + 7r~7r~) < 7.4 x 10~ 7 
and B {B- -> K-R--K+) < 4.2 x 10~ 7 . 

In conclusion, we present searches for the standard 
model suppressed B meson decays B~ — ► K + Tt~Tt~ and 
B~ — > K~ K~ir + . We do not see any evidence of these 
decays and obtain improved 90% confidence level upper 
limits on the branching fractions. These results super- 
sede those of our previous publication [To| and can be 
used to constrain models of new physics. 
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ing organizations that support BaBAR. The collaborat- 
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FOM (The Netherlands), NFR (Norway), MES (Russia), 
MEC (Spain), and STFC (United Kingdom). Individuals 
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FIG. 3: Construction of the confidence region bands. The left (right) plot is for B~ — + K + tt~it~ (B~ —* K~ K~-k + ). In each 
figure the blue dotted line shows the expected central value of Nat as a function of the true branching fraction, the green solid 
(red dashed) lines show the 90% confidence level upper and lower limits including statistical and systematic errors (statistical 
errors only), the black dashed horizontal line marks the position of the previous upper limit [ljj, and the black dash-dotted 
lines indicate the results of this study. 
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